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Abstract 

Several lattice collaborations performing simulations with 2+1 light dynamical quarks have 
experienced difficulties in fitting their data with standard Nf = 3 chiral expansions at next- 
to- leading order, yielding low values of the quark condensate and/or the decay constant 
in the Nf = 3 chiral limit. A reordering of these expansions seems required to analyse 
these data in a consistent way. We discuss such a reordering, known as Resummed Chiral 
Perturbation Theory, in the case of pseudoscalar masses and decay constants, pion and 
kaon electromagnetic form factors and form factors. We show that it provides a good 
fit of the recent results of two lattice collaborations (PACS-CS and RBC/UKQCD). We 
describe the emerging picture for the pattern of chiral symmetry breaking, marked by a 
strong dependence of the observables on the strange quark mass and thus a significant 
difference between chiral symmetry breaking in the Nf = 2 and Nf = 3 chiral limits. We 
discuss the consequences for the ratio of decay constants Fx/F n and the form factor 
at vanishing momentum transfer. 



1 Introduction 



Recent improvements in lattice simulations have paved the way for unquenched simulations with 
three light dynamical (or sea) quarks. Even though these simulations are still at an early stage, their 
results are already very instructive, especially when one is interested in the low-energy dynamics of 
strong interactions. Indeed, there is a deep connection between these first-principle simulations and 
the effective description of the light hadronic degrees of freedom, i.e., Chiral Perturbation Theory 
(xPT) [H [2] . The ability for the lattice to investigate the quark-mass dependence of hadronic observ- 
ables opens the possibility to determine some poorly known low-energy constants. Conversely, 
the chiral expansions derived within xPT can be used to extrapolate or interpolate lattice data down 
to the physical light quark masses. 

It turns out that several collaborations performing simulations with 2+1 dynamical quarks re- 
ported difficulties when fitting their data with Nf = 3 next-to-leading-order (NLO) chiral expansions 
for pseudoscalar masses, decay constants, and -fQ 3 form factors [21 El E] , forcing some of the collab- 
orations to rely on Nf = 2 chiral expansions only (for instance using on Heavy-Kaon %PT to perform 
their extrapolations [7]). Moreover, some chiral order parameters (in particular the quark condensate 
and the decay constant) were seen to decrease significantly when one moved from the Nf = 2 chiral 
limit where m u , rrid — > but m s is kept at its physical value [1] to the Nf = 3 chiral limit where 
ra u ^m^m s — > [2j. Defining the chiral order parameters in these two chiral limits (denoted lirnpfy 
for simplicity): 

S(JV» = -\im Nf (0\uu\0) , F 2 (N f ) = \im Nf F% , (1) 
PACS-CS [3] and MILC ^ quote for instance: 

PACS-CS : = 1.089 ± 0.045 , = 1.245 ± 0.010 , (2) 

F{3) £(3) 

MILC : M = 1 . us±a06 (^), H = 1JB±M r(^). 

Assuming that there are no intrinsic problems with the lattice data, these results can be interpreted as 
the fact that the Nf = 3 chiral series do not converge quickly, because their leading-order term in the 
chiral counting is not numerically dominant and competes with (formally) higher-order contributions. 
Calculations up to NNLO are clearly useful to settle the issue [101 [TT], [121 E3] - An exploratory study 
of the RBC/UKQCD data [H [5] on form factors and Fx/F n (the ratio of the kaon and pion 
decay constants) showed that NNLO terms are not negligible though of the expected size, the NLO 
terms of some observables being in this particular study anomalously small for some known reasons. 
However, this study did not consider the RBC/UKQCD result for the decay constant F n itself, which 
does not exhibit the quark-mass dependence predicted by chiral perturbation theory [HE]. 

Quite interestingly, a similar pattern with suppressed Nf = 3 chiral order parameters seems to 
emerge from the experimental data on irn and irK scattering, indicating a good convergence of Nf = 2 
chiral series for pion observables [151 El El El El [201 ED E2] , but difficulties in Nf = 3 chiral series, 
even once NNLO terms are included and modeled through resonance saturation [23[ I24j. As recalled 
in refs. [10[ [25l I26| . the situation seems to be well under control in the tttt system, but the fits yield 
sometimes contradictions once irK data is considered. It was advocated in refs. [271 [28] that such a 
situation could indeed occur in the low-energy dynamics of QCD and it was highlighted there that 
the strange quark may play a very special role, due to its light mass of order 0(Aqcd)- Significant 
vacuum fluctuations of ss pairs may lead to instabilities in the convergence of expansions, where 
instabilities are defined as a numerical competition between the terms considered as leading and next- 
to-leading in the chiral counting. This effect would be related to a large violation of the Zweig rule 
in the scalar sector, indicated by values of the 0(p 4 ) LECs L4 and Lq significantly different from 
specific ("critical") values [291 [30] ( see a ^ so re f s - [SUES] for reviews of the subject and a discussion 
of the baryon sector). 
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A pessimistic way of considering the problem would consist in dismissing the whole xPT as soon 
as problems of convergence arise. A less restrictive point of view was adopted in ref. [33], assuming 
that: 

• only some ("good") observables have convergent expansions, when expressed in terms of the 
couplings arising in the chiral Lagrangian. 

• a series is considered as convergent when the sum of LO and NLO terms is large compared to 
the remaining part of the series. 

• the resulting formulae must be treated analytically, without neglecting higher-order corrections 
when reexpressing low-energy constants in terms of observables. 

This framework was coined Resummed Chiral Perturbation Theory (Re^PT) for reasons that will 
be recalled later in the present article. In particular, a specific prescription for the unitarity pieces, 
described in sec. 12.21 will be added to this set of assumptions, defining precisely how we choose to 
split the chiral series of the chosen observables between leading, next-to-leading and higher orders in 
our framework. 

The recent 2 + 1 dynamical simulations provide new and relevant information on the impact 
of ss fluctuations related to the presence of strange quarks in the sea. Conversely, Re%PT can 
provide a more appropriate treatment of chiral extrapolations if the hints of suppressed Nf = 3 
quark condensate and decay constants are confirmed. In principle, this analysis would require lattice 
data performed with several u, d, s quark masses (whose renormalized values are known) and transfer 
momenta (in the case of form factors, scattering amplitudes), but where the continuum and infinite- 
volume limit have already been performed (a — > 0,L — > oo). Unfortunately, such data sets are 
not (yet) available. Some collaborations [e.g., MILC [8j [9]] provide numbers directly in the physical 
limit, performing the chiral extrapolation at the same time as the continuum limit. This prevents one 
from testing different alternatives concerning chiral extrapolation, even though the results sometimes 
contradict the usual assumptions of xPT (for instance concerning the size of the quark condensate 
and the decay constant). Others [BMW [331 EE! , ETMC [36], TWQCD-JLQCD [37]] do not provide 
the decay constants and the renormalized quark masses mandatory for such a study. Finally, some 
collaborations [RBC/UKQCD || |], PACS-CS [3]] have performed their analysis only at one 
particular lattice spacing and/or one particular volume, without estimating the systematics associated 
with the continuum and infinite-volume limits fully. 

In view of this situation, we will restrict ourselves to the studies made by RBC/UKQCD and 
PACS-CS. The fact that only statistical errors are quoted in both cases prevents us from using a 
full-fledged statistical treatment [25], but we hope that our limited study will provide some incentive 
for lattice collaborations to apply the same framework to their data when performing chiral extrapo- 
lations. In Sec. [2] we provide the basics of our procedure with the illustration of the electromagnetic 
form factor of the pion, and we recall the results obtained for the masses and decay constants. In 
Sees. [3] and IU we apply the same formalism to the electromagnetic form factor of the kaon and the 
Kis form factors. In Sec. [U we consider the same observables in the case of lattice simulations for 
non-physical values of the quark masses. In Sec. [6l we fit our chiral expansions to different sets of lat- 
tice data, and analyse the emerging pattern of Nf = 3 chiral symmetry breaking, before concluding. 
An appendix summarises the main features of the results obtained by the two lattice collaborations 
which are analysed in this article. 
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2 General arguments and the electromagnetic form factor of the 
pion 

2.1 Expansion of "good" observables 

As explained in greater details in refs. [381 E31 EH] , if one expects a numerical competition between LO 
and NLO chiral series of quantities of interest, one cannot perform expansions on arbitrary functions 
of such quantities. Let us assume for instance that an observable A has the following chiral expansion: 

A = A LO + A NLO + Ad A (4) 

The statement of a good convergence of the series correspond to 5 A much smaller than 1, but it does 
not involve necessarily that Alo is dominant numerically with respect to the rest of the series. 
One has the formal chiral expansion for any observable F = f(A): 

F = F LO + F NLO + F5F , (5) 
Flo = f{A LO ) F NLO = A NL of'{A L o) , (6) 

SF = 1--^ —— [A -Au>-MA]. (7) 

Depending on the value of Xa = Alo /A, i- e - 5 the saturation of the chiral expansion of A by its first 
term, the chiral series of F may or may not converge well. One has in particular the two limiting 
behaviours: 

X A ^l SF^-ll^AdA, (8) 

y _>n ap^i f'^A^f'^AXA (q\ 

Xa ^° 6F ^ 1 -m~JiA) A+ JiA) A6A - (9) 

In the first case, a bound on 5 A implies a bound on 5F, meaning that F converges well provided 
that A does. But in the second case, the size of 5F is by no means driven by that of 5 A. As an 
illustration, if we take the "observables" B = 1/A and C = VA, we obtain: 

SB = (10) 
5C = 1 ^--JXa" 1=5 A . (11) 

Even if we set 5 A = 0, we would need Xa = Alo I A above 41% to ensure that \5B\ < 10%, and Xa 
above 76% to ensure that |<5C| < 10%. 

Therefore, if the chiral expansions of some of the observables considered are not saturated by 
their leading-order term, we cannot take an arbitrary function of these observables, consider its 
chiral expansion, and assume that it will converge. In such a situation, one has therefore to select 
the right set of observables for which converging series can be written. Actually, this statement 
of convergence is equivalent to state that the generating functional of QCD is well reproduced at 
low energies by xPT. The quantities that are assumed to have convergent chiral series from the 
starting point of chiral perturbation theory are the QCD correlators of axial and vector currents 
as well as pseudoscalar and scalar densities. In addition, the convergence is expected to be good 
only away from the singularities (poles, cuts . . . ) corresponding to resonances and channel openings. 
Observables involving pseudoscalar mesons as external states will be obtained by applying the LSZ 
reduction formula on correlators involving axial currents, which will yields additional factors of the 
pseudoscalar decay constants. One can easily state some of the observables associated to each 0(p 4 ) 
LEC that are simple to determine on lattice simulations and/or to extract from experiment: 
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. L X ,L 2 ,L 3 : K u form factors (F^F K f,F^F K g,F^F K h), 

• L4,L^: tt,K decay constants (F^,F^), 

• L 6 ,L 8 : 7T, K masses (i^Af* F| M|), 

• L7: 7/ mass and decay constant (F^,F^M^), 

• Lg\ Pion electromagnetic form factor (F^Fy), 

• L\q: it — > ewy form factor (F^Avap) and r spectral functions (VV — AA). 

We will define the expansion of a "good" observable through the following procedure: 

• We take the chiral expansion of the observable in terms of the couplings of the effective La- 
grangian (low-energy constants, or LECs: Bo,Fo,L>i. . . ). 

• We replace the pseudoscalar masses at leading order by the physical ones only when physical 
arguments indicate that the convergence of the series will be improved. In practice we perform 
this replacement so that the nonanalytic structures imposed by unitarity are located at the 
physical poles, thresholds. . . 

• We keep track of the remainders (collecting NNLO contributions and higher orders) explicitly 
and treat the resulting chiral expansions as algebraic identities, without performing further 
expansions. 

Following this procedure (and taking the isospin limit m u = = m), observables can be 
expressed in terms of LO quantities: 

m-^. m = ^, 

as well as NLO LECs and remainders. The first two quantities in eq. (|12p are of particular relevance, 
since they express two main order parameters of Nj = 3 chiral symmetry breaking, the quark 
condensate and the pseudoscalar decay constant, in physical units. They also assess the saturation of 
the chiral expansion of F%M% and F% by their leading order. The third quantity measures the relative 
size of the quark masses in a framework where the strange quark is supposed to play a peculiar role 
in the chiral structure of QCD vacuum. In the following, we will also use the quantity ^(3): 

y(S) = *2* = *W (13, 
V ' Ml Z(3) V 7 

which assesses the saturation of the chiral expansion of the pion mass by its leading order. Up to now, 
we have considered the masses and decay constants of Goldstone bosons [281 EQj EH] and observables 
derived from tttt scattering [381 [33] in this framework. 

Exploiting the fact that some quantities describing the dynamics of pseudoscalar mesons are well 
measured and inverting the relationships between these observables and LECs, we can express the 
0(p 4 ) LECs in terms of 

• Masses, decay constants, form factors. . . 

• The three leading-order parameters described in eq. (112p 

• The remainders associated to each observable, assumed to be small (convergence). 
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These expressions can be exploited in the chiral expansions of other "good" observables, in order to 
express the latter quantities in terms of LO quantities and remainders only. The comparison with 
experimental information should then provide more information on the pattern of chiral symmetry 
breaking in the Nj = 3 chiral limit. 

As shown in refs. [38j ESI EE], this procedure applied to masses and decay constants allows one 
to resum higher-order contributions in chiral series from L4 and Lq low-energy constants, which 
encode the effect of ss pairs on the structure of the chiral vacuum. It may induce a significant m s - 
dependence in the pattern of chiral symmetry breaking and can generate a numerical competition 
between LO and NLO in Nf = 3 chiral series. This feature is related to the name of Resummed 
Chiral Perturbation Theory chosen to describe this particular treatment of chiral expansions. This 
framework is compatible with the usual treatment of chiral series in the limit where the latter are 
saturated by their LO term, but it allows for a consistent treatment of the series even if there is a 
significant competition of LO and NLO contributions for some of the observables. 



2.2 Electromagnetic form factor 

We will illustrate our procedure with the example of the electromagnetic pion form factor 

(7r+|i M |7r+) = (p + p')^y(*), 



(14) 



where the electromagnetic current is j M = Vjf + V^f/\/3, p (p') is the momentum of the incoming 
(outgoing) pion, and t = (p' — p) 2 . As explained in Sec. 12.11 we obtain this form factor in xPT from 
((A^ + )^ A^ + j^) , leading to the product F 2 Fy through the LSZ reduction formula (the pion decay 
constant F n stems from the wave- function renormalisation) . In the case of the electromagnetic form 
factor, "good" observables are thus obtained from F 2 Fy at low energies away from singularities, i.e., 
the right-hand cuts, starting from t > AM 2 . 

We obtain the following bare expansion, in agreement with Refs. [39 1 1401 [41] : 



F 2 Z(3) + M 2 Y(3)[8(r + 1)L\ + 8L r 5 
1 



+t 



2LI 



32vr 2 

, 2 



M 2 Y(3) 



41og^f + (r + l)lo£ 



o 2 

M K 
M 2 



o 2 

1, M n I, M K 1 
- log + - log — f + - 

3 /X z fl D 



+i[t - 4M 2 Y(3)]J n7T (t) + l[ t _ 2(r + l)M 2 Y(3)]J KK (t) + RF 2 F v (t) , 



(15) 



with the nonanalytic pieces from the two-meson channels encoded in the J function [2] . Mp denotes 
the leading order of the chiral expansion of 0(p 2 ): 



o 2 



M^= M^y(3) 



M K = M 2 J-±^Y{3) 



M v =M 2 ±(2r + l)Y(3) 



(16) 



and Y(3) is the ratio defined in eq.(|13|). We have added §iF 2 Fy(t), a polynomial function of t 
collecting remainders: 



^F 2 Fy(t) = (MF 2 Fy) + —(^Fy), + 0(t 2 ) , 



(17) 



with {^F 2 F V ) = 0(m 2 ) and ($lF%Ffi)i = 0(m q ). 

At next-to-leading order in the chiral expansion, the nonanalytic dependence on quark masses 
and momenta arises through the unitarity function Jpq. Following our prescription, we compute 
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the functions J (and J = J — sJ'(O)) with the physical values of M 2 , Mj^, M 2 , rather than their 
leading-order expansion, i.e. we define the chiral expansion in RexPT as eq. (|15p with: 



Jpp{t) 



16tH 



dx 



4M 2 X (X 



4M| 



1 



16vr 2 



crlog 



a 



1 



CJ + l 



+ 2 



4M 2 

i-V 1 - < 18 > 



so that our expansion of the form factor eq. (|15p features a unitarity cut from the two-pion channel 
starting at t = AM 2 (id. for the two-kaon channel). Indeed, from general arguments of unitarity, we 

o 2 

know that the higher-order corrections will shift the start of the right-hand cut from 4 Mp to AM p. 
Unfortunately, unitarity does not provide us more information on the structure of the cut (and in 
particular the coefficient multiplying the J function) due to the perturbative nature of the chiral 
expansion. 

When needed, we will obtain the RexPT expansion of other observables by performing the same 
replacement for the functions J and J occurring in the definition of the loop integrals Kpq , L pq and 
M P Q in ref. [2] . However, we do not perform any further replacement neither in the unitary functions 
nor in the rest of the expressions: for instance, we have not modified the functions multiplying the 
J functions, nor the chiral logarithms coming from the tadpole terms in eq. (I15p . since we have no 
way of determining if the latter modifications would improve or spoil the convergence of the series 0. 
Adopting a dispersive point of view, we can say that the position of the cuts are imposed by unitarity, 
but not the value of the induced imaginary parts at low energies and that of the subtraction constants 
(polynomials). 

One checks easily that the NLO chiral expansion of the electromagnetic form factor in ref. [39] 
can be recovered: 

F${t) = 1 + 2H^(t) + H KK (t) + 0(p 4 ) , (19) 



with 



Hpp(t) = —2 



-(« 

12 v 



AMp) J PP (t) - - 



t 1 



6 32vr 2 



log 



Ml 
M 2 



+ 1 + 



288tH 



2t 
3K 



jr 



(20) 



In the case where Fo is small compared to F n , as hinted at by lattice simulations and NNLO fits of 
chiral series |26tl42|. F%Fy(t) is expected to exhibit a better convergence than Fy(t) in our framework 
according to Eqs. (I19p -(|20p. Similar expressions hold for other observables: good observables will 
generally come multiplied by powers of physical pseudoscalar decay constants (one for each external 
pseudoscalar meson involved). 



2.3 Pion electromagnetic square radius 

The electromagnetic square radius of the pion is the low-energy observable associated with Fy\ 

,dFi 



QFt 



dt 



■(0) 



(21) 



Following the previous discussion of the form factor, the product of F 2 and (r 2 )y is the quantity 
expected to exhibit a good convergence in our framework, eq. (|15p yields the corresponding expansion 
of (r 2 ) v : 



(r 2 Yv 



6 

J2 



2ALc 



32vr 2 



I + 1y(3) + 
6 9 w 18M 



M 2 

72-^ + 1)^(3) 



K 



+ (r 2 T v e v 



(22) 



1 This procedure is slightly different from the approach taken in Ref. [25], where this substitution was performed 
everywhere in the unitarity functions J, K, L, M and in the tadpole logarithms. It turns out that the difference is 
usually very small: the unitarity functions yield only a small contribution below the first threshold, and there is only 
a logarithmic difference in the case of the tadpole. 
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where we have introduced the scale-independent combination ALg = L r ^{ji) — L g (/i), with: 



32^ 



1. M« 
6 \i A 



+ 



12 



log 



o 2 
M_K 



(23) 



and 



6 Qggg 



F 4 



(24) 



The pion electromagnetic square radius is well-determined, and it is expected to suffer only mildly 
from higher-order corrections being an observable involving pions. We will thus use this observable 
to express Lg: 



F. 



AL 9 = ^-(r 2 r v [l-e* v ] + 



1 



12 



32vr 2 



1 1 

12 + 9 y(3) + 36M 



M 2 

-^-(r + l)y(3) 



A 



(25) 



Equivalent relations for other LECs, namely £4,5,6,8 will be discussed later in eqs. (|39l) - (H2l . Lg(M p ) 
can thus be estimated as a function of r, Y (3) and NNLO remainder. For instance, if we take Y(3) = 1 
and r = 2M^/M 2 — 1 (corresponding to LO estimates holding in the case of a fast convergence) and 
the central experimental value 



{r 2 )\r = 0.451 ± 0.031 fm 2 , 



(26) 



we obtain! L r Q (M p ) = 6.77-10" 3 in the ball park of usual estimates of this LEC, such as (6.9±0.7)-10^ 3 
at 0(p 4 ) in Ref. [2], (5.93 ± 0.43) • 10" 3 at <3(p 6 ) in Ref. @T]. On the other hand, if Y(3) tends to 
zero while the radius and its NNLO remainder remain finite, we see that Lg becomes very large (we 
will see later that the other electromagnetic radii remain also finite in this limit). 



2.4 Masses and decay constants 

Exactly as in the previous section, we can write the identities for Fp and FpMp in terms of 
r,X(3),Z(3), the 0(p 4 ) LECs L4, L5, Lq, L§ and remainders: 



FlZ{3) + MlY{3) [8(r + 2)U A + 8L r 5 ) 

M 2 Y(3) 41og^ + (r + l)log^ 



1 

32^ 



+ FX , 



F 2 K = FlZ(3) + MlY{3) [8(r + 2)L r A + 4(r + 1)L\ 

o 2 o 2 



(28) 



(29) 



1 



32vr- 



■MlY(3) 



F 2 M 2 = F£M*X{3) + M*[Y{3)Y [16(r + 2)L£ + 16L£ 

o 2 o 2 

1 



(30) 



32vr- 



:M 4 [y(3)]^ 



o 2 

,g^ + (r + l)log^ + i(2r + l)log^ 



2 In this article, we take the following numerical values for the pseudoscalar masses, the pion decay constant and the 
renormalisation scale: 



Mk = 0.13957 Gev , M K = 0.4957 Gev , M„ = 0.5478 Gev , = 0.0922 Gev , /x = 0.770 Gev . (27) 
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= -(r + 1) F 2 M 2 X(3) + M*[F(3)] 2 [16(r + 2)U e + 8(r + 



(31) 



32vr 



-M*[Y(3)] 5 



o 2 o 2 o 2 

- log + -(r + 1) log — f + — (2r + 1) log -J- 



A* 



+ F 2 K M 2 K d K ■ 



The pion and kaon masses are well known experimentally. As far as the decay constants are 
concerned, F n and Fk are accessible at a high precision through leptonic decays (7172 and Ki 2 respec- 
tively [33]) which provide in the framework of the Standard Model (SM) 



\V US /V ud \ x F K /F n = 0.2758 ± 0.0005 , 



(32) 



which can be combined with the very accurate determination of the first element of the CKM matrix 
V u d from super-allowed + — > + nuclear beta decays [35] 

(33) 
(34) 



V ud = 0.97425 ± 0.00022 , 
the unitarity of the CKM matrix and the smallness of the \ V U \,\ matrix element to get: 



F-r 



92.2 ± 0.3MeV 



\SM 



1.192 ±0.006, 



\SM ~ « Fk/F^I 

In the following we take the value of F n in eq. (|34p . expecting the deviation from this SM determi- 
nation to be rather small. On the other hand, we will not fix the value of Fx/F n , keeping it as a 
free parameter of the fit to the lattice data. A deviation from its value eq, (|34p would hint at physics 
beyond the Standard Model contributing to flavour-changing charged currents, in addition to the 
usual V — A term from the W bosons. 

Similarly to the case of Lg and the electromagnetic square radius of the pion and following |33[ [25] 
we will then invert the relationships eqs. (|281) -(|3ip in order to reexpress the four NLO LECs in terms 
of X(3), Z(3),r, the pion and kaon masses, the pion decay constant and the ratio Fk/F^: 

1 Fl 



y(3)AL 4 
Y(3)AL 5 



8(r + 2) M% 



[1 _ n (r) _ Z (3) - e] 



(35) 
(36) 



y 2 (3)AL 6 



y 2 (3)AL 8 



F 2 



16(r + 2) Ml 
1 F 2 



[l-e(r)-X(3)-d\ 



16 MV 



(37) 
(38) 



ALj = LKfi) — Li(fi) combine the (renormalized and quark-mass independent) constants 1/4,5,6,8 & n d 
chiral logarithms so that they are independent of the renormalisation scale \x: 

o 2 



327r 2 L 4 (/i) 



32tt 2 L 5 (^) 



32vr 2 L 6 (^) 



32vr 2 L 8 (/i) 



7: log 



o 2 

M K 



1 



8 n 2 8(r-l)(r + 2) 



o 2 

(4r + 1) log + (2r + 1) log — f 



log 



o 2 
M_K 

M 2 



+ 21og^ 



+ 



1 



8(r - 1) 



1 

16 
1 

16 



o 2 



o 2 

log — f + - log 

o 2 o 2 

. M K 2, Af„ 
log — f + - log 



M K M n 

o 2 



M K 



16 (r + 2)(r-l) 



M K 



+ 



16(r - 1) 



o 2 



o 2 

„, M K . M 
3 log— ^ + log 

O ^ 



'/ 



o 2 



(39) 



(40) 



(41) 



(42) 
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The four equalities eqs. (|35H38]) are exact, since they are a mere reexpression of the bare chiral 
series for F%, F^, F%M% and F^M^. d,d' and e,e' are combinations of remainders associated with 
the chiral expansions of tt, K masses and decay constants respectively: 



d = r -^-d n - ( e(r) + -?—) d K , d' = d-d T , (43) 

r — 1 y r — 1 J 

r + 1 / 2 \ 

e = -e^ - r)(r) H ) ex , e' = e - e n . (44) 

r — 1 \ r — 1 y 

d', e' are quantities of 0(mm s ) whereas d, e scale like O(m^). In addition, the right hand-side of 
these equations involves the r-dependent functions: 

The properties of these equations, and in particular, the fact that they lead to a resummation of 
L4 and Lq contributions once inserted in the chiral expansion of other observables, were discussed at 
length in refs. [25J ESI [38] . Note that contrary to the case of Lg, these LECs always appear multiplied 
by powers of V(3) in chiral series, since they correspond to operators with one or two powers of the 
scalar source in the chiral Lagrangian, and arise in chiral expansions always multiplied by one or 
two powers of Bq. Therefore, when Y(3) tends to zero, Li will exhibit a logarithmically divergent 
behaviour (due to the tadpole logarithms) which will not affect observables though. 



3 Kaon electromagnetic form factors 

The method described in the previous section can easily be generalized to other observables. Of 
particular interest are the kaon electromagnetic form factors and the Kir form factor which will be 
discussed in the following sections. 

3.1 Definition 

The kaon vector form factors [39^ |4"T] are defined as: 

(K+\j IM \K+) = (p + pYF{f + (t), (K^K ) = (p+pYFfit), (46) 

with the same convention as in the case of the pion electromagnetic form factor. All of them are 
associated with the P-wave projection of the crossed channel. Following the discussion in O we 
expect F^Fy + and F^-Fy° to have good convergence properties away from the singularities (opening 
thresholds. . . ). Expanding these form factors and reexpressing some couplings in terms of r,Y(3) 
and Z(3), we obtain the bare expansion of the vector form factors: 

2 ° 2 
-^Ff{t) = - I 92^1og^f- I 2^[t-4M^(3)]J 7r7r (i) (47) 

+ -^[t - 2(r + l)M*Y(3)]J KK {t) + ^F 2 K Ff(t) , 
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7?2 

-f2 F V (*) 



Z(3) + ||[F(3)] [8(r + 2)L^ + 4(r + 



32vr 2 F 2 



F(3) 



o 2 o 2 

^og^ + ^(r + l)log^f + i(2r + l)log^ 



t 



IV, 



32vr 2 



/r 2 

o 2 



o 2 



1 . M n 1 , M x , 1 

^ lo e — + o lo g —f- + g 



6 ° ^ 2 



+^2 [* - 4M 2 y(3)] J_(t) + J^[t - 2(r + l)M 2 Y(3)]J KK (t) 
+^F 2 K F^ + (t), 

where tftF^Fy (t) and ^Si.F^Fy + (t) are polynomial functions of t collecting remainders: 



(48) 



KFlFf{t) 



{KFlFf) 1 + 0(e 



F 2 



RF. 



P K 



?2 t?K+ 



2 i?K+^ 



(49) 
(50) 



with (?R.Fy)o = 0(m 2 ) and (MFy)i = 0{m q ). We have divided the expressions of the form fac- 
tors eqs. (I48| )- (147p by a numerical factor F 2 for sole purpose of convenience, in order to deal with 
dimensionless quantities. 

In the limit where all the observables are saturated by their leading order, the standard NLO 
chiral expansions of the vector form factors |391 W\\ can be recovered by expanding the ratio F 2 / F^ 
at next-to-leading order and replacing the leading order masses by the physical ones: 

Ff(t) = -H^t) + H KK {t) + 0(p 4 ) , F# + (t) = F v (t) + Ff(t) + 0{p A ) , (51) 

with Hpq defined as: 



H PQ (t) 



1 

F? 



1 
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t — 2Spq + 



A 2 v 



JpQ{t) -^r JpQ{t) -l kpQ + 2^ 



+ ^1, (52) 



involving Epq = Mp + Mq and Apq = Mp — Mq, and Jp Q {t) = Jp Q {t) - tJ' PQ (t). 
3.2 Kaon electromagnetic radii 

In a similar way to the pion form factor, the K + electromagnetic square radius is given by 



6 

K 



2ALc 



1 "'^^^.^H.™ 

M„ K 



32vr 2 6 



+ (r^fef , (53) 



with the remainder: 



K+ _ 6 (KF 2 F^H 
e v 



2 >r 



(54) 



Replacing ALg by its value in terms of the pion radius, eq. ()25|) . leads to the following relation: 

o 2 

? 2 i2\K+ f-t „K+< 



F^)ra-er)-F^rTv(l-e VJ 



1 



32^ 



log^f + ^ 

M.,r 



Ml 



3M| 



(r + l)y(3) 



(55) 
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where the right-hand side is a very small correction for any reasonable value of r and Y{3), so 
that the electromagnetic square radius of the charged kaon is essentially predicted to be (r 2 )y + — 



F'£/F'k x (r 2 )y — 0.32 fm 2 . Two experiments have measured this radius, leading to the average |43j : 



(r 2 )f = 

The square radius of the neutral kaon reads: 

1 



0.314 ± 0.035 fm 2 



with the remainder: 



32vr 2 



log 



o 2 

Mr 

o 2 

Mtt 



+ - 3 Y(3) 



Ml 
3M 2 



(r + l)y(3) 



P K° 
e V 



6 (KF£ )i 



F 2 ( r 2) 



2\AT° 
V 



The current experimental average is 

(r 2 )£° = -0.077 ± 0.010 fm 2 . 
Eqs. (|55p and (|57|) yield the following relation between the electromagnetic radii: 



F 2 
K 

Fl 



(<r 2 >ni-en-<r 2 >f(i-ef)) . 



(56) 



(57) 



(58) 



(59) 



(60) 



which is fulfilled using the experimental values of the radii and the SM value of Fx / F n (the remainders 
must be on the large side of their allowed value according to the dimensional estimation discussed in 
Sec. I5.3|) . In principle, the knowledge of these remainders and the determination of the radii with a 
high precision would allow us to determine Fk/F w accurately. 



4 Kit form factors 
4.1 Definition 

Among the quantities that can be determined from lattice simulations, one can single out the 
form factors defined as: 

V2(K + \u^s\7r°) = (p'+ P rf + (t)+(p'- P ru(t). (6i) 

/+ corresponds to P-wave projection of the transition, whereas its 5-wave comes from 

/„(*) = / + (t) + -^/_(t), (62) 

where Apq = Mp — Mq. Following the discussion in El F- K Fxf+ and F n Fxf- are expected to have 
good convergence properties away from the singularities (opening thresholds. . . ). Exactly as before, 
their chiral expansions can be expressed in terms of r,X(3),Z(3), NLO low-energy constants (L4, 
L5 and Lg) and remainders. Reexpressing L4 and L5 using eqs. ([35|) -([36 ]) yields the following bare 
expansions of the form factors: 

F v F K f + (t) = F " ± ^ + \[tM^(t) + tM r Kv {t) - L Kw (t) - L Kv (t)} (63) 
+2tU 9 + F n F K d + + te + , 
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F„F K f_(t) 



t?i _ TP? 



-{Ml - Ml)[M^{t) + M r Kv (t)} 
5(t - Ml - M 2 K ) + |(r + 3)Mjy(3) 



(64) 



3(t - Ml - Ml) + -(r + 3)A4 2 y(3) 

-2(A4 - M^)Lg + F^fk (d_ - d+) + t(e- - e+) , 

where d± = 0(m 2 ) and e± = 0(m q ) combine the remainders from the form factors and the decay 
constants: 



(MF n F K f + ) 



F^F K {d„-d + ) = (KF n F K f„) + 

F n F K e + = (XEkFkU)!, 
F 7T F K (e.-e + ) = (ftF n F K f + ) u 



Ffa + F 2 K e K 
2 

Ffa - F 2 K e K 



(65) 

(66) 

(67) 
(68) 



where the remainder F n Fx(^if±)(t) are defined as before. When performing our fits to lattice data, 
we will also express Lg in terms of the pion radius using eq. (|25]). Inserting eqs. (|63p - (|64j ) into eq. (|62j) 
leads to the following expression for the scalar form factor: 



F*F K f (t) 



Fl + Fl 



+ 



2 
t 



~ ^ ~ ~ 2 L KTr{t) ~ -L Kr) (t) 



t 



Artt 

-K K „(t) 



(69) 



4A 



■K Kv {t) 



5(t - Ml - Ml) + -(r + 3)MlY(3) 
3(t - Ml - Ml) + i(r + 3)M*Y(3) 



+(F 7T F K d+ + te+) ( 1 - + (F n F K d_ + te * 



Artv 



In the limit where all expansions are saturated by their leading-order contribution, the well-known 
expression for the vector form factor is recovered: 



(70) 



as well as that for f^ n [39]. 



4.2 Callan-Treiman point, its soft kaon analog and the form factor at zero mo- 
mentum transfer 



According to the Callan-Treiman theorem [36], in the soft-pion limit (p' 2 = Ml = 0), the scalar form 
factor at t = A Kn = Ml - Ml (Callan-Tr eiman point) should be equal to Fx/F n . This implies that 
FjcFirfoiA-Kn) ~ Ffc van i snes i n the Nf = 2 chiral limit m u = rrid = m — > 0. There is a soft-kaon 
analog of this theorem holding at t = Ax n = —Ak-k, stating that FxF n fo(AKn) — F% vanishes in 
the Nf = 3 chiral limit. At these particular points, eq. (169p reads: 



F n F K f (A K « 
1 



Fl 



+ -^Kktt{^-Ktt) 

-jK Kv (A Kit ) 



- -Lktt(Aktt) - -Lkt){Ak-k) 

-10Ml + ^-(r + 3)MlY{3) 



(71) 



-5Ml + -(r + 3)MlY(3) 



+ F n F K d_ + A Kn e_ , 
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F n F K f (-A KjT ) = Fl 

— jKktt(—^Ktt) 



-10M 2 K + -(r + 3)M 2 Y(3) 



-6M 2 K + -(r + 3)M 2 Y(3) + F n F K (2d, 



(72) 



A Kn (2e + - e_) 



One can check explicitely that these expressions fulfill the Callan-Treiman theorem and its soft- 
kaon analog (the K and L contributions canceling each other) provided the following constraints on 
the NNLO remainders 

d- = 0(mm s ) , e_ = 0(m) , (73) 

meaning that d_ and e_ are suppressed compared to d + and e+. 

We can define the discrepancies from the Callan-Treiman theorem(s): 



A CT = f (A Kn ) - ^ 



A CT = fo(-A K n) - ■ 



(74) 



These NLO quantities can be expressed from eqs. ()7ip -(|72 p . embedding the fact that Act is 1/r- 
suppressed compared to Act- For comparison, these quantities have been calculated in standard 
XPT at one-loop order in the isospin limit [47] : 



A 



CT 



-3.5 • 10" 



Act = 0.03 . 



(75) 



It has in fact been shown in refs. \48\ |4"9] that a precise assessment of the scalar form factor at the 
Callan-Treiman points could probe physics beyond the Standard Model in the strange quark sector, in 
particular right-handed couplings of quarks to W bosons. The pioneering work [IS] led to a reanalysis 
of Kps data by several collaborations [5Q |. 151 } 152]. which at present show a good/marginal agreement 
with the Standard Model except for the NA48 collaboration [50] exhibiting a 4.5a deviation still 
unsolved. 

The Kits vector form factor at zero momentum transfer is another quantity of interest. Indeed, the 
measurement of decays can be analysed in the framework of the Standard Model to determine 
the product |V^ S /+(0)|, and thus the CKM matrix element \V US \. A recent fit to \V U( i\ (from super- 
allowed + — > + nuclear decays), |V^ s |/+(0) (from K^), and \V us /V U( i\Fk / F^ (from -k^ and Kgz) 
together with the unitarity of the CKM matrix led to 



/+(0) 



SM 



0.959 ± 0.005 , (76) 
with a strong correlation between these two 



and a value of Fk / F w | SM in full agreement with eq. ([31 
quantities. 

Deviation of /+(0) from this value would be an indication of new physics, so that this quantity 
plays a particularly important role to test the Standard Model in the light quark sector. A direct 
determination of these quantities on the lattice as well as a well-controlled method to extrapolate 
lattice data down to the physical quark masses are naturally crucial to get a proper assessment of 
the uncertainties (from statistical, but also systematic origins). 



5 Observables for lattice simulations at different quark masses 

As explained in the previous section, we can use the relations eqs. (f35|) - (f36|) (decay constants), 
eqs. (|3Tl) - (j3"8l) (masses), eq. (f25l) (pion electromagnetic square radius)... to express NLO LECs in 
terms of r, X(3), Z(3), accurately measured observables and remainders. These relations can be 
inserted in the chiral expansions of other observables (such as kaon or form factors, or meson- 
meson scattering), which can be used to constrain r, X{3) and Z(3). For each new observable, one or 
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several remainders are introduced, which are assumed to be small but nevertheless limit the accuracy 
of the chiral series. 

As indicated in the introduction, one can also consider lattice simulations, where the same observ- 
ables are considered at different values of the quark masses. The interest is twofold. First, the lattice 
simulations probe the m s -sensitivity of observables, hard to estimate from continuum measurements, 
but with deep connection with the pattern of Nf = 3 chiral symmetry breaking. Second, lattice 
simulations [RBC/UKQCD 01 13 El, PACS-CS [3]] have encountered difficulties in their fits of NLO 
Nf = 3 chiral expansions. Let us remark that the inclusion of NNLO terms for form factors 
and Fk/F-x to fits of the RBC/UKQCD data seems to solve convergence issues for these particular 
quantities assuming no Zweig-rule violation in the scalar sector [42]: a good \ 2 is obtained with a 
rather good convergence of the chiral series with NNLO terms of the expected size (this work did 
not discuss F n itself, for which the problems of convergence seem the most acute [4[ [14] ) . 

Fitting these data will offer us the opportunity to extract relevant information on chiral symmetry 
breaking and to check the consistency of our picture concerning the numerical competition between 
LO and NLO terms. We consider simulations with 3 dynamical flavours (m, rh, m s ) and denote X 
the values for the lattice quantities (and X the corresponding value for physical quark masses). We 
introduce the ratios: 

m s m 
P = — . 1 = — > ( 77 ) 



m. 



in* 



in addition to the ratio of physical quark masses r and the chiral parameters arising in the leading- 
order Lagrangian in eqs. ([T2l) and (Ti~3j) . 



5.1 Masses and decay constants 

Proceeding as before in this new setting, we obtain the following expansions for the decay constants: 



F 2 M 2 

-f 2 = Z{3) + -j£pqrY(3) 



Ml 



1 



Fl 32tt 2 



pqrY(3) 



4 log 



- + 2 



M 2 



+ I - + 1 ] log 



2 

Mk 
M 2 



+ 



(78) 
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(79) 



F 2 



where the LO contributions to the simulated pseudoscalar masses are involved: 

-2 ~2 /i \ ~ 2 



M w = pqrMlY{3) , 
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+ l MlY(3) 
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v 



+ l MlY(3) , (80) 



and ep are remainders of 0(rh 2 ) [fh q denotes either rh s or rh). We have divided by the physical value 
of F 2 in order to deal with dimensionless quantities. In a similar way, we obtain the bare expansions 
of the masses: 
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j^=P<ir\x(3) + 1 £pqrlY(3)] 2 
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(82) 



where dp are remainders of 0(m 2 ). We have divided by the physical value of F 2 M 2 in order to 
deal with dimensionless quantities. As explained before, we use eqs. (f35|) - ([38|) to express the (mass- 
independent) chiral couplings £4,5,6,8 m terms of r,X(3),Z(3) and the physical masses and decay 
constants. 



5.2 K e3 form factors 

We obtain for the lattice vector form factor: 

Fj K f + {t) ■- 



H±Ik + hfM^{t) + tM r Kv (t) - L Kn (t) - LkM 
+2tL r 9 + F n F K d + + te+ , 



and the scalar form factor: 
F v F K f (t) ~ 
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(83) 



(84) 



where Lpq, Kpq, Mpq are evaluated with the leading-order pseudoscalar masses at the simulated 
quark masses using eq. (|80p . apart from the Jpq function which is evaluated at the simulated ("phys- 
ical") pion and kaon masses using eqs. (|8ip -(i82l). In the above formulae, the decay constants on the 
right-hand side arise from the reexpression of L4 and L5, and should be understood as a short-hand 
notation of the full expressions in eqs. (|78p - (179| ). For the vector form factor, we can trade Lg for the 
pion electromagnetic square radius using eq. ([251) . 



5.3 Remainders 

The expressions for the simulated masses, decay constants, form factors and electromagnetic square 
radius involve unknown remainders. These remainders collect all the contributions coming from 
NNLO, NNNLO and higher orders. They can be evaluated by resonance saturation [53], involving a 
hadronic scale Ah only mildly affected by the actual value of the quark masses (mass of the p, K* . . .). 
In order to keep track of the scaling of the remainders with the quark masses, we take the following 
NNLO estimates which involves the hadronic scale at the fourth power: 



d,e,d K ,e K ,d+ = O 



d' e' d_ = O 



( 2M%M 2 K 
d n = d— d' 



O 



O 



F 2 M 2 
A 4 

A 4 



O 



6 



M 2 K 



(r 2 r v A% 



(85) 

(86) 
(87) 



15 



where M 2 and M|- follow the known dependence of the remainders on m and m s , whereas F 2 is 
inserted when a dimensionful constant with no dependence on m q is required. 

We can use the known scaling of the remainders to perform their extrapolation to the simulated 
quark masses: 

d n = p 2 d — p 2 qrd' , 

2 2/ 

e n = p e — p qre , 
d+ = p 2 d+ , e+ = pe + , 

6 Fit to lattice values 

6.1 Data and parameters 

We are now in a position to build the \ 2 function to be minimised for the sets of data that we will 
consider. The inputs are the following ones, as recalled in App. [Aj 

• Fl,F^FlMl,F^M\ for known values of the quark masses (m,m s ) (RBC/UKQCD and 
PACS-CS) 

• F n Fxf+ and F n Fxfo for several transfer momenta (RBC/UKQCD) 

• We consider the quantities given by these collaborations corresponding to light quark masses 
and small momenta ("Subset" fit) where chiral perturbation theory is valid. 

In Appendix B we will also consider all the quantities available even though some points at 
higher masses and momenta might certainly be outside the region of validity of Re%PT ("Total" fit) 
to check the stability of our results as well as to illustrate the interest of having more data points 
to determine the remainders more accurately. Obviously, we would need more points for lower pion 
and kaon masses. 

The uncertainties on these quantities were obtained by combining the uncertainties in quadrature: 
no correlation between the various observables is provided in the articles of both collaborations, and 
we have only the statistical errors (no estimate of the systematic uncertainties was available for the 
quantities of interest here). The parameters entering the fit are: 

• Quantities from the leading-order chiral Lagrangian X(3), Z(3), r 

• NNLO remainders: d, d' ', e, e' (in all cases), d + , e+, cL, e_, e\ (for a fit to decay constants), 

• The value of the ratio of decay constants Fx/F n , 

• The value of p re f = vn S)Te j jm s for a lattice set of reference, providing the equivalence between 
lattice and physical quark masses (when possible). 

The last quantity is estimated by both collaborations, but we found it interesting to keep this param- 
eter free in the fit, in order to take partially into account systematic effects related to lattice spacing. 
Since the quark masses are expressed in a mass-independent scheme involving only multiplicative 
renormalisation, we can determine the value of p = rh s /m a for any lattice set once we know p for a 
given reference set using m s = (rh s /rh Sjre f) ■ p re f- We had to fix the last two parameters in the case 
of RBC/UKQCD "Subset" fits, due to the limited number of unitarity (unquenched) points available 
for the masses and decay constants (only 2 different pairs of quark masses). 

When computing the values of the observables from NLO chiral expansions, we need the values 
of the masses and decay constants for the simulated quark masses (for instance in the unitarity 



V F,M n J y 2 V 2 
e x =(^)V( e -!^Ve<), (89) 
d- = p 2 qrd- , e_ = pqre- . (90) 
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functions J). In such a case, we computed systematically the values of the decay constants and 
masses from their chiral expansions (|78p -(|82p. rather than plugging in their "measured" values on 
the lattice. This distinction may have some importance for the form factors eqs. (|63p and (I69p . 
where we have reexpressed L4 and L5 in terms of F% and F^, but where the latter quantities stand 
for their chiral expansion in terms of LECs (the fits to RBC/UKQCD would be slightly improved 
compared to the ones presented here if we used the measured values of F% and F\ rather than the 
computed ones). 

In addition, the mass and decay constant F„ of the 77 are needed for the evaluation of the 
loop integral Jpq (and related unitarity functions). They are obtained at a sufficient accuracy for 
such purposes using the two following LO formulae reminiscent of the Gell-Mann-Okubo relation: 

F" = pK-\^, F^ = -^FlMl-\FlMl. (91) 

We constrain the remainders in the ranges indicated in Table [5j Once the (MINUIT-powered) 
fit has converged, we can estimate a large body of quantities: NLO LECs, Nf = 2 chiral order 
parameters, values of the scalar form factor at zero momentum transfer, at the Callan-Treiman 
point and its soft kaon analog, test of the convergence of the series. We have propagated the er- 
rors exploiting the covariance matrix provided by MINOS, assuming that all uncertainties follow a 
Gaussian distribution. 

6.2 Discussion 

Our results are summarised in Table [TJ The first series of rows corresponds to the outcome of 
the fit, whereas the lower rows are quantities derived from the results of the fit (LO LECs, NLO 
LECs, quantities in the Nf = 2 chiral limit, quantities, relative fraction of LO/NLO/remainders 
contributions at the minimum for several observables) , and the last row is the \ 2 P er degree of 
freedom. 

The "Subset" fit of RBC/UKQCD results includes only a limited number of data points for 
the masses and decay constants, which forces us to fix one of the parameters of the fit, namely the 
simulated strange quark mass El and to impose some bounds on the size of the higher-order remainders, 
based on a simple estimate from resonance saturation described in App. |Bj Indeed, some of these 
remainders are pushed to the limits of their range when the set of data is too small, because there 
is not enough information for MINUIT to choose a particular value for these remainders (keeping 
them free would lead to a larger contribution from higher orders and to a further decrease of the 
leading-order contribution). The situation improves when higher masses and momenta are included 
(see the fits to "Total" sets presented in App. [C"|) . so that the remainders remain small, within the 
window discussed in App. |Bj It is interesting to notice that the results are consistent with the fits 
to PACS-CS data ("Subset" or "Total"): r is close to 25 (i.e close to the ratio 2M^/M% - 1 even 
though its value was left free in our framework), the quark condensate remains around X(3) ~ 0.5, 
the squared decay constant is Z(S) ~ 0.6, leading to a value for the LO term of the squared pion 
mass F(3) ~ 0.8. 

The ratio of decay constants Fx/F^ (left free in our fit) comes out slightly larger (smaller) for 
PACS-CS (RBC/UKQCD) than its Standard Model value eq. (|34p in the fits of the masses and decay 
constants. We obtain also values of simulated strange quark masses and of the physical strange 
quark mass in good agreement with the results obtained by the two collaborations (the discrepancy 
between RBC/UKQCD and PACS-CS is due to the different choice of renormalisation procedure, 
which explains the low value obtained by the PACS-CS collaboration |54|). 



3 Letting all parameters free gives comparable results for the central values, but some parameters get very large 
uncertainties, larger than their allowed range. Propagating the errors in such a situation would be meaningless, and 
reporting the results of this fit would not provide much more information than the constrained fit that we present here. 
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PACS - CS Subset 
Without 


RBC/UKQCD Subset 
With Kta 


r 

X(3) 

y(3) 

Z(3) 

Rem. at limit 

rh s ,ref/m s 


26.5 ±2.3 
0.59 ±0.21 
0.90 ±0.22 
0.66 ± 0.09 
1.237 ±0.025 
none 
1.24 ±0.08 


23.2 ± 1.5 
0.20 ±0.14 
0.43 ±0.30 
0.46 ±0.04 
1.148 ±0.015 
d, e 
1.15* 


m s {2 GeV)[MeV] 
m(2 GeV)[MeV] 
S (2 GeV)[GeV] 
FofMeV] 


70 ±4 
2.6 ±0.3 
3.34 ± 1.18 
74.8 ± 4.9 


107 
4.6 ±0.3 
0.92 ±0.67 
62.2 ±2.5 


l 4 (m) • io 3 
l 5 (m) • io 3 

L 6 ( M ) • 10 3 
L B (ji) ■ 10 3 
£g(/x) • 10 3 


-0.1 ±0.2 
1.8 ±0.4 
0.1 ±0.4 
0.8 ±0.7 

X 


2.4 ±2.0 
1.8 ± 1.6 
4.7 ± 7.1 
4.4 ± 7.1 
4.4 ±2.8 


X(2) 
Y(2) 
Z{2) 
B{2 GeV)[GeV] 
F[MeV] 

h 
h 

B/B 
F/F 


0.90 ±0.01 
1.04 ±0.02 
0.87 ±0.02 
3.83 ± 0.50 
85.8 ±0.7 
5.0 ±2.1 
4.5 ±0.5 
1.51 ±0.51 
1.15 ± 0.26 
1.15 ±0.08 


0.90 ±0.02 
1.00 ±0.03 
0.90 ±0.02 
2.09 ±0.19 
87.7 ±0.8 
-0.6 + 3.7 
3.3 ±0.5 
4.52 ±2.83 
2.28 ± 1.39 
1.41 ±0.06 


/+(0) 
A C r • 10 3 
A' CT ■ 10 3 


1.004 ± 0.149 

X 
X 


0.985 ± 0.008 
-0.2 ± 12.1 
-126 ± 104 


(rX° N 2 ] 


X 
X 


0.248 ± 0.156 
-0.027 ±0.106 


F n F K f+(0) 


0.66 + 0.22 + 0.12 
0.44 + 0.48 + 0.08 
0.60 + 0.30 + 0.10 
0.42 + 0.50 + 0.08 

X 


0.45 + 0.69 - 0.14 
0.34 + 0.76 - 0.10 
0.20 + 0.95 -0.15 
0.14 + 0.97-0.11 
0.40 + 0.75 - 0.15 


X 2 /N 


0.9/3 


4.4/8 



Table 1: Results of fits performed on the data from the PACS-CS [3] and RBC/UKQCD [H [6] 
collaborations on pseudoscalar masses and decay constants, and form factors in the case of 
RBC/UKQCD. In all cases, we considered only data with light pions (Subset) and only statistical 
errors are shown. In the RBC/UKQCD case, we fixed the lattice strange quark mass (marked with 
a star). The LECs are given at the scale [i = m p . In the PACS-CS case, the form factor at zero 
momentum transfer is a prediction of the fit (with an error combining those obtained from the fit and 
the maximal contribution allowed for the remainder from dimensional estimation). The penultimate 
set of rows collects the relative fractions of LO/NLO/remainders for decay constants, masses and 
K13 form factor at vanishing transfer momentum (for RBC/UKQCD) at the minimum. 
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The decay constant in the Nf = 3 chiral limit is found to be rather low, in agreement with other 
recent works \26\ W2\ . The pattern of Nf = 3 chiral symmetry breaking (with low quark condensate 
and decay constant) is reflected by the values obtained for the low-energy constants L4 and Lq, 
which are both positive and do not show any sign of Zweig suppression. Such large values of L4 
and/or Lq have been obtained in several earlier works: dispersive analysis of scalar form factors |29j . 
dispersive treatment of ttK scattering [23], J ftp decay into a vector meson and two pseudoscalars |55j 
(with a value of Lq compatible with zero), preliminary NNLO Nf = 3 fits to pseudoscalar masses, 
decay constants, Ku decay and irK scattering data [26]. . . Large values of L4 and Lq are known to 
induce a significant dependence of the chiral order parameter on the strange quark mass, and it is 
not surprising to witness a strong suppression from the Nf = 2 chiral limit to the Nf = 3 one. Let 
us notice that we obtain values for Nf = 2 chiral order parameters @ which are in agreement with 
the hypothesis of standard confirming that this framework is indeed appropriate in the Nf = 2 
sector as shown by the recent data on decays [201 [21], though this does not seem to be the case 
for Nf = 3. The other LECs £5,8,9 have values in agreement with conventional estimates (this was 
expected in particular for Lg since our framework induces modifications that are only sub-leading for 
vector quantities such as the pion electromagnetic form factor). 

Let us notice some specificities of our treatment of the chiral expansions. Tadpoles diagrams 
generate chiral logarithms of the form Mp log (Mp/fj, ) which can prove quite troublesome to fit. For 
instance, the RBC/UKQCD collaboration [U \5\ [6] finds a better agreement of their data on decay 
constants with polynomial fits than with chiral series. In our treatment, these chiral logarithms always 

o 2 

involve the leading-order mass Mp- Therefore, the limit of a small Y(3) tames the chiral logarithms 
in our expansions, so that these logarithms become hard to distinguish from a polynomial at the 
numerical level on the range of masses where xPT could be valid. Furthermore large contributions 
from NLO LECs, and in particular from L4 and Lq as just discussed, will enhance the quadratic 
dependence on the quark masses and thus our chiral expressions will mimic a polynomial dependence 
on the quark masses that cannot be reproduced in the more usual treatment of chiral expansions. 
These mechanisms could explain why chiral logarithms are often difficult to identify in lattice data, 
in addition to other effects (heavy strange quark mass, lattice systematics. . . ). 

The K13 form factor at zero momentum transfer, /+(0) = /o(0), involves only LECs related to 
decay constants and masses, eq. f|63j) . In principle, it can be predicted from a fit of the latter quantities 
up to the determination of the remainder c£ + . We quote the corresponding results in Table [Tj where 
the central value for /+(0) corresponds to remainders set to zero. The uncertainty on this quantity 
includes the maximal size allowed for the remainders d+ based on dimensional estimate (see Table E]) , 
as well as the uncertainties coming from the parameters of the fit. Clearly, the NNLO remainder 
cL|_ hinders any accurate determination of /+(0), unless their value is also precisely determined from 
the fit, which is possible once data on form factors themselves is included (Tabled]). The values 
obtained for /o(0) are somewhat larger than the Standard Model value eq. (|76j) . as well as those 
obtained from the RBC/UKQCD collaboration using different forms for the extrapolation in quark 
masses [5] [6] . This illustrates the importance of the mass extrapolation for lattice simulations at the 
level of accuracy aimed at currently. A particular attention was paid in ref. [6] to the structure of 
the chiral expansion of /+(0) = 1 + fi + f^, where /2 is the NLO contribution, which involves only a 
combination of chiral logarithms divided by Fq\ 

h(x) = 1 + - 2X „ log x . 
1 — x L 

(92) 

fi is often said to be free from LECs and thus known precisely from Chiral Perturbation Theory. 

4 The expressions for the Nf = 2 order parameters X{2), Y(2), Z{2) can be obtained by setting rh = 0,m s = m s in 
the expressions for F% and F%M%, eqs. (|78[) and (|81[) . The matching expressions for £3,4 in terms of £4,5,6,8 differ from 
the usual ones 2, 56 by factors of Y(2) and Y(3) that are easy to recover. 



h 



256tt 2 F 2 



{M 2 K + Ml)h 



Ml 

Ml 



+ {Ml + M 2 v )h\ 



Ml 
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This statement is not totally correct for the following reasons. One usually assumes that the value 
of Fq is close to that of F n , so that it can be replaced in actual calculations by the physical value 
of the pion decay constant leading to the estimate f2 — —0.023. Since the difference between the 
two quantities is a higher-order effect, one can always perform this replacement. However, one has 
to determine how large the NNLO term can be with such a prescription, and consequently how 
well the chiral series for / + converges. If Fq is significantly lower than F n as indicated not only by 
our fit, but by other recent estimates [261112], the convergence is expected to be rather slow, forcing 
us to treat the NLO contribution to /+(0) more carefully. We advocated that correlators of vector 
and axial currents yields observables with good convergence properties, selecting F^F^foiO)- In this 
case, we should replace Fq by F w Fk in the evaluation of eq. (192ft . as can be checked in our expression 
for /+, eq. @3}. 

Once the form factors are included in our fits, Lg can be determined even though the fit does 
not constrain for this particular LEC tightly. In Table [TJ the deviations from the Callan-Treiman 
relation at t = Ax n and its soft kaon analog at —A^n are given. Their values are of the expected 
size for SU(Nf) chiral-symmetry breaking quantities for Nf = 2,3 flavours respectively, and thus 
compatible with the one obtained in standard xPT. The values of the square radii of the charged 
and neutral kaons, also shown have rather large uncertainties and are thus within the experimental 
error bars. 

In the last lines of our tables, we have indicated for each fit the contribution from LO, NLO 
and remainders to pseudoscalar decay constants and masses for values of the parameters at the 
minimum of the fit. We can see that the series converge well on overall (remainder much smaller 
than LO+NLO), but that the LO term is far from saturating the series. The values of 3^(3) obtained 
is smaller than 1, reducing the contribution from chiral logarithms compared to that from the NLO 
LECs. We can compare these results with those from a fit of the same observables, where the 
NLO and higher contributions (chiral logarithms /ip, LECs Li, remainders) are computed replacing 
2mBo, (to + m s )BQ and Fq by the physical pion and kaon masses and the pion decay constant. This 
is exactly equivalent to performing the same fit as before with the following replacements in the NLO 
and higher-order contributions: 



M 



K 



Ml 



1. 



Ml 



2M 2 K 



Mr 



p 



2M 2 K 



Ml 



2M 2 K 



rj(r) ->■ 77(7-0), e(r) ->■ e(r ) - 2X(3)- 



ro 



1' 



log 



Mr 

o 2 

M P 



Y(3) 



log 



Ml 
M 2 : 



1. 



(93) 



(94) 



both for the observables that we consider, eqs. (|78p -(|84 p . and the equations allowing the determi- 
nation of £4,5,6,8,9, eqs. (J25J and (|59 > (132" ]) . For PACS-CS, this leads to x 2 /N = 1.1/3 (compared 
to our result 0.9/3), with very similar values for the fundamental parameters r,X(3),Y(3),Z(3). 
For RBC/UKQCD, the fitting procedure yields % 2 /N = 9.5/8 (compared to our result 4.4/8), with 
much more uncertain values of the fundamental parameters (e.g., r = 14.9 + 12.1, X{3) = 0.30 + 0.26, 
Y(3) = 0.68+0.60). This is not particularly surprising since our fits to the PACS-CS data led to values 
of r and Y(3) in good agreement with eq. (|93p . but not the RBC/UKQCD ones. The corresponding 
convergence of the pseudoscalar masses and decay constants is then (the relative contribution for 
LO, NLO and higher orders is given here): 



PACS - CS [NLO phys. masses] 
RBC/UKQCD [NLO phys. masses] 



F 2 

F 2 
K 



F 2 
K 



0.64 + 0.26 + 0.10, 
0.42 + 0.51 + 0.07 
0.45 + 0.70 - 0.15 
0.34 + 0.76 -0.10, 



FlMl 
F 2 K M 2 K 
FlMl 



0.67 + 0.24 + 0.09, 
0.50 + 0.44 + 0.06, 
0.31 + 0.81 - 0.12, 
0.15 + 0.94-0.11. 



There is no saturation of the series by their leading order. We see that our formulae yields results 
that are in good agreement with those obtained after reexpressing the NLO contributions in terms 
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of F n , M w ,Mx in the PACS-CS case, where 1^(3) is close to 1. On the other hand, when Y(3) is not 
close to 1 (for instance in the RBC/UKQCD case), our formulae provide more efficient and accurate 
fits (lower x 2 , smaller error bars). From a more methodological point of view, we avoid a perturbative 
reexpression of LECs in terms of F n , M n , Mk in a regime where it is not justified. 

These trends can be compared interestingly with the fits done by the lattice collaborations them- 
selves, with a different treatment of the chiral series than ours. For instance, the MILC collab- 
oration [9] observed from fits with staggered chiral perturbation theory that M% received NNLO 
corrections of the same size as NLO contributions, canceling each other to a large extent, with small 
NNNLO corrections (the latter being taken as analytic in quark masses and lattice spacings), whereas 
F n exhibited no problems of convergence. On the other hand, the RBC/UKQCD collaboration [4] 
experienced difficulties in fitting F n both in Nf = 2 and Nf = 3 theories. They also noticed that 
fits to Mj£ and Fk using the Nf = 3 chiral expansion led to very significant NLO contributions (of 
order 50%) when data up to the kaon mass scale was included, and they conclude that higher-order 
corrections could be very significant (up to 30%). 

At this point, we should emphasize that our framework does not contain any bias concerning the 
size of X(3), V(3) and Z{3) or on the relative size of the LO and NLO contributions. It is compatible 
with the usual assumptions that chiral series of decay constants, squared masses. . . are saturated by 
their LO contribution, but it can also accommodate situations where there is a numerical competition 
between LO and NLO terms. It turns out that the lattice data set from the RBC/UKQCD and 
PACS-CS collaborations favour values for the three quantities X{3), Y(3), Z(3) smaller than 1, with 
a x 2 /d.o.f. which ranges from fairly good to excellent. Our results confirm the difficulties reported 
by the two collaborations to fit Nf = 3 NLO chiral expressions, and highlights the improvement 
provided by our Re%PT formulae for the extrapolations in quark masses of these quantities. As a 
further check, we have performed fits where we have taken that the physical masses (and not the LO 
ones) in the unitarity functions J, K, L, M and the argument of the chiral logarithms (similarly to 
what was done in ref. [25]). The quality and parameters of the fits are almost unchanged, and the 
outcome for the derived quantities is also very similar, meaning that the relevant issue is the proper 
choice of the "good observables" whose chiral series converge well. 

7 Conclusion 

Recent lattice simulations with light 2+1 dynamical fermions have encountered difficulties to fit their 
results for pseudoscalar masses, decay constants and form factors with chiral expansions obtained 
from Nf = 3 Chiral Perturbation Theory at next-to-leading order. Such fits of poor quality can 
be related to the fact that chiral series are not saturated by their leading order, so that there is 
a numerical competition between leading-order contributions - from the decay constant and/or the 
condensate in the Nf = 3 chiral limit (m u = = m s = 0) - and next-to-leading-order contributions 
- in particular from L4 and Lg, related to the Zweig-rule violation in the scalar sector, enhanced by 
m s and not accurately known. 

If there is such a competition, one must decide which observables are expected to have a good 
overall convergence (small higher-order contributions). According to the assumed equivalence of the 
%PT and QCD generating functionals at low energies, it seems reasonable to consider observables 
derived from correlators of axial and vector currents as well as pseudoscalar and scalar densities, 
as done here. For these observables, one must treat chiral series with a particular care, avoiding 
the perturbative reexpression of LECs in terms of observables while neglecting higher orders (this 
can be easily done by introducing remainders corresponding to NNLO and higher contributions) and 
choosing how unitarity contributions should be treated to define the structure of the chiral expansion 
and its splitting into leading, next-to-leading and higher-order terms. Such a set of prescriptions was 
introduced some time ago under the name of Resummed Chiral Perturbation Theory. In the present 
paper, we have recalled the basic ingredients of this framework and applied it to observables related 
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to pseudoscalar masses, decay constants, and kaon and pion form factors (electromagnetic and 
ones). This allowed us to illustrate how 0(p 4 ) LECs £4,5,6, 8, 9 can be reexpressed in terms of the 
leading-order quantities X(3), Z(3),r as well as experimental values of observables (pion and kaon 
decay constants and masses, square electromagnetic radius of the pion) and associated remainders. 

Then we have turned to 2+1 lattice simulations where these observables were obtained for several 
sets of quark masses: PACS-CS (decay constants and masses only) [3] and RBC/UKQCD (decay 
constants, masses, form factors) [HG2 EJ- We performed fits to data corresponding only to light 
quark masses and small momenta, but checked the stability of our procedure by considering also fits 
to all data available (unitary points). Since only statistical uncertainties (without correlations) are 
publicly available for each of the points, we performed naive fits with Gaussian errors, in order to 
determine the leading-order parameters of the chiral Lagrangian as well as NNLO remainders and 
the ratio of decay constants. 

The fits are generally of a good quality, with a good consistency when one compares subsets 
coming from the same collaboration. This allows one to determine the values of the LO quantities 
as well as the NNLO remainders, with a good accuracy in the case of PACS-CS, with a more limited 
precision for RBC/UKQCD because of the restricted number of low-mass points. One observes that: 

• The decay constant and the quark condensate in the Nf = 3 limit (m u = = m s = 0) are 
both small and suppressed compared to the Nf = 2 case {jn u = md = and m s physical). 

• The low-energy constants L4 and Lq do not follow the Zweig rule suppression generally advo- 
cated to set them to zero. 

• The other low-energy constants L5, L$ and Lg have values in good agreement with previous 
estimates. 

• The ratio of quark masses r remains quite close to the most simple estimate from pseudoscalar 
masses. 

• Nf = 2 chiral order parameters are in good agreement with the values extracted from 
decays. 

• When the sets of data are large enough, the NNLO remainders remain in the expected range 
from a naive dimensional estimate. 

• The expected numerical competition between LO and NLO chiral expansions indeed occurs for 
Fl F*, FlMl and F*M*. 

Beyond this description of the pattern of Nf = 3 chiral symmetry breaking and its implication 
for the convergence of chiral expansions, we can also make a few predictions. The values obtained 
for the kaon electromagnetic radii are in good agreement with experimental data. In the case of 
RBC/UKQCD, the value obtained for /+(0) with our fits is slightly larger than the ones quoted by 
the collaboration, relying on alternative formulae for the chiral expansion of the form factors. 
This has naturally an impact on the determination of | V us | , considering the level of accuracy achieved 
in decays [Hj . 

Lattice simulations are able to investigate the dependence of observables on quark masses, which 
makes them very valuable tools to investigate the chiral structure of QCD vacuum. Conversely, 
any improvement in our understanding of chiral symmetry breaking will help reducing systematics 
associated with chiral extrapolations in lattice determinations. In light of the discussion presented in 
this paper, it would be very helpful that more lattice collaborations present the dependence of their 
results on the quark masses and/or study alternative ways of performing the extrapolation down 
to physical quark masses in order to assess the related systematics precisely. It would also be very 
interesting to have lattice results for observables related to the scalar channel, and thus difficult to 
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determine from experiment. For instance, the pion and kaon scalar form factors, simulated on the 
lattice and analysed in the framework of Resummed Chiral Perturbation, would provide an interesting 
complement to the present discussion. 
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A Lattice inputs 

A.l RBC/UKQCD Collaboration 

We first consider the RBC/UKQCD Collaboration simulations with 2+1 dynamical flavours [H[5l|6] 
performed with domain-wall fermions at one lattice spacing a^ 1 = 1.729(28) GeV. The calculations 
are performed on two volumes, 16 3 x 32 and 24 3 x 64 ((2.74) 3 fm 3 ) (with a fifth dimension of length 
16), at each quark mass, except the lightest mass which is only simulated on the larger volume. 
They performed a non-perturbative renormalisation to relate the lattice quark masses to those in the 
RI-MOM scheme. 

The only points that we will use are those where the sea and valence quark masses are identical. 
There are four sets corresponding to such a situation for pseudoscalar masses and decay constants 
in Ref. [1], corresponding to a(m lat — m res ) and a(m^ a * — m res ) being respectively (0.005,0.040), 
(0.010,0.040), (0.020,0.020), (0.030,0.030). where am res = 0.00315(2). The quark masses are given 
in the RI-MOM scheme, but they can be related to the MS scheme through a multiplicative factor 
m(2 GeV) = Z m a~ 1 (am lat ) that drops in all the input quantities (which involve only ratio of quark 
masses in the same setting). 

We obtain the following values expressed in units of 10~ 3 GeV -2 for the pseudoscalar masses 
and decay constants [UET]. The uncertainties here are purely statistical and do not include those 
induced by the uncertainty on the value of the lattice spacing. 



Masses 


(p, q) 


F 2 


F 2 M 2 


Fl 


F 2 K Ml 


Subset 


(0.005,0.040) 
(0.010,0.040) 
(0.020,0.020) 
(0.030,0.030) 


(1.15,0.189) 
(1.15,0.304) 
(0.616,1) 
(0.883, 1) 


10.98 ±0.16 
12.85 ±0.16 
15.58 ± 0.34 
17.82 ± 0.36 


1.196 ±0.022 
2.249 ± 0.036 
4.851 ±0.107 
8.038 ±0.166 


14.11 ±0.19 
15.59 ±0.18 
15.58 ±0.34 
17.82 ±0.36 


4.644 ± 0.076 
5.730 ± 0.082 
4.851 ±0.107 
8.038 ±0.166 





In two papers [5j[6], the RBC/UKQCD collaboration investigated the form factors /o and 
/+ using in particular twisted boundary conditions to obtain a sample transfer momenta, with the 
same two sets of values corresponding to nondegenerate masses (a(rh lat — m res ),a(m ! s a! — m res )) = 
(0.005, 0.040), (0.010, 0.040). 

The set with the lighter u, d quark masses yields the following values: 



t 


60.7 


59.87 


38.1 


21.6 


0.30 


F*F K f (t) 


12.68 ±0.17 


12.73 ±0.17 


12.49 ±0.17 


12.32 ±0.17 


12.15 ±0.16 


F w F K f+(t) 


X 


X 


12.71 ±0.176 


12.42 ±0.175 


12.15 ±0.17 


Subset 






★ 


★ 
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r 


28.8 ± 0.4 ± 1.6 


rh, rP f/m q 


1.150 


Fk/Fk 


1.205 ± 0.018 ± 0.062 


m s (2 GeV)[MeV] 


107.3 ±4.4 ±9.7 ±4.9 


m(2 GeV)[MeV] 


3.72 ±0.16 ±0.33 ±0.18 


5(2 GeV)[GeV] 


2.52 ±0.11 ±0.23 ±0.12 


F[MeV] 


81.2 ±2.9 ±5.7 




3.13 ±0.33 ±0.24 


4 


4.43 ± 0.14 ± 0.77 



Table 2: Results obtained by the RBC/UKQCD collaboration in ref. @]. 



t 


-44.00 


-129.3 


-204.9 


-389.2 


FnF K fo(t) 


11.68 ±0.21 


10.95 ± 0.32 


10.77 ±0.23 


9.667 ±0.28 


F«F K f + (t) 


X 


X 


X 


X 


Subset 











The set with the heavier value leads to the following values for the scalar form factor: 



t 


35.42 


-90.51 


-195.3 


-205.0 


-385.2 


F^F K f (t) 
Subset 


14.28 ±0.17 


13.05 ±0.21 


11.64 ±0.38 


12.89 ± 0.40 


11.83 ±0.75 



We have considered also a subset of data, indicated with stars, where the convergence of chiral 
series is expected to be particularly good. This amount to considering only non-degenerate u, d, s 
quark masses, and to drop the points for t < —0.2 Gev -2 [i.e., the two points for /q in ref. [5j 
corresponding to the lowest values of transfer momentum]. 

The values of the physical quark masses (m and m s ) and the lattice spacing are obtained by 
studying the dependence of the mass of tt, K and Q, hadrons on these three parameters and tuning 
them to reproduce the physical hadron masses. If we call rh s ^ re f the value of the strange quark mass 
corresponding to the set (0.005,0.040), the RBC/UKQCD collaboration obtained rh Sjre f/m s = 1.150. 
Considering the uncertainty associated with such a determination (in particular the role played by 
the form of the chiral extrapolation used for n and K), we will not assume this value in our fit, but 
rather include this quantity as a parameter of our fit, and scale the other ratios involving a simulated 
strange quark mass over the physical value. 

Fits to the Nf = 2 and Nf = 3 NLO chiral series for pseudoscalar masses and decay constants 
were performed in ref. [I]. It turned out that the Nf = 3 chiral expansions led to rather poor fits 
(large x 2 P er d.o.f), in particular for decay constants, unless they put stringent cuts on the values of 
quark masses where such expansions should hold. This led the authors in ref. [3] to perform fits to 
Nf = 2 NLO chiral expansions. In addition, in ref. [14] . NNLO SU(2) chiral expansions were shown 
to have only a limited utility to extrapolate the data: many more data points would be needed to 
fix the size of the combinations of 0(p 6 ) counterterms involved The results obtained in ref. [4] that 
are relevant for our discussion are summarised in table [2 

In addition, two different values for /+(0) were obtained in refs. [51 [6] from the same gauge 
configurations, using either data for the scalar form factor or data for both form factors, and applying 
a pole ansatz based on either Nf = 3 or Nf = 2 chiral perturbation theory for form factors [7]: 

/+(0) = 0.964 ± 0.033 ± 0.0034 ± 0.0014 |5] , /+(0) = 0.960(11) © . (95) 

A. 2 PACS-CS collaboration 

The PACS-CS collaboration [3] has investigated the pseudoscalar masses and decay constants with a 
large sample of light quark masses, for one particular value of lattice spacing a" 1 = 2.176(31) GeV, on 
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a 32 3 x 64 volume. They used a non-perturbatively 0(a)-improved Wilson quark action and performed 
the renormalisation of quark masses perturbatively at one loop (with tadpole improvement), with 
the following results: 



( arn ud \a,mf s ) 


(p,q) 


F 2 

x TT 


FX 


F 2 
K 


F\Ml 


Subset 


(0.001,0.040) 
(0.006, 0.041) 
(0.010,0.036) 
(0.011,0.042) 
(0.021,0.045) 
(0.031,0.047) 


(1.410,0.040) 
(1.456,0.138) 
(1.256,0.271) 
(1.519,0.260) 
(1.577,0.466) 
(1.663,0.652) 


10.19 ± 1.09 
11.51 ±0.26 
12.48 ±0.21 
13.25 ±0.18 
17.23 ±0.73 
19.09 ±0.51 


0.247 ±0.035 
1.007 ±0.031 
1.846 ±0.041 
2.242 ± 0.036 
5.595 ±0.239 
9.397 ±0.254 


14.29 ±0.48 
15.49 ±0.22 
15.37 ±0.16 
16.83 ±0.23 
19.87 ±0.65 
21.01 ±0.58 


4.385 ±0.151 
5.459 ± 0.088 
5.200 ±0.067 
6.791 ± 0.096 
10.11 ±0.33 
13.08 ±0.37 





Once again, the uncertainties are of purely statistical origin, and they do not include the uncertainty 
coming from the determination of the lattice spacing. We have also considered a subset of data, 
indicated with stars, where the convergence of the chiral series is expected to be particularly good. 
This amounts to considering the three lightest values of the pion masses. 



r 


28.8 ±0.4 


y(3) 


0.88 ±0.01 


Z(3) 


0.76 ±0.04 


F K /F n 


1.189 ±0.020 


rh S! ref/m s 


1.19 


m s (2 GeV)[MeV] 


72.72 ± 0.78 


m(2 GeV)[MeV] 


2.527 ±0.047 


B {2 GeV)[GeV] 


3.869 ± 0.092 


F [MeV] 


83.8 ± 6.4 


l 4 (m) • io 3 


-0.06 ±0.10 


lm ■ io 3 


1.45 ±0.07 


L 6 (p) • 10 3 


0.03 ±0.05 


L 8 ( M ) • io 3 


0.61 ±0.04 


Y(2) 


0.96 ±0.01 


Z(2) 


0.88 ±0.01 


B(2 GeV)[GeV] 


0.96 ±0.01 


F[MeV] 


88.2 ± 3.4 


h 


3.14 ±0.23 


h 


4.04 ±0.19 




1.205 ±0.014 


B/B 


1.073 ± 0.055 


F/F 


1.065 ± 0.058 



Table 3: Results obtained by the CP-PACS collaboration with one-loop perturbative renormalisation 
and extrapolation to the physical limit [3]. The values for the quantities in the Nf = 2 chiral limit 
correspond to Nf = 2 fits to the so-called Range I with finite-size effects included. 

The values of the physical quark masses (m and m s ) and the lattice spacing are obtained by 
studying the dependence of the mass of tt, K and O hadrons on these three parameters and tuning 
them to reproduce the physical hadron masses. If we call rh s ^ re f the value of the strange quark mass 
corresponding to the set (0.0016,0.0399), we obtain m Sjre f/m s = 1.19. Considering the uncertainty 
associated with such a determination (in particular the role played by the form of the chiral extrap- 
olation used for tt and K), we will not assume this value in our fit, but rather include this quantity 
as a parameter of our fit, and scale the other ratios involving a simulated strange quark mass over 
the physical value. 
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Fits to the Nf = 2 and Nf = 3 NLO chiral series for pseudoscalar masses and decay constants 
were performed in ref. [3]. It turned out that the Nt = 3 chiral expansions led to rather poor fits, 
related to very significant NLO contributions compared to LO terms, in particular for the decay 
constants, related to large contributions from kaon loops. In other words, the dependence of these 
quantities on the strange quark mass seen in these simulations is not accounted for properly by NLO 
SU(3) chiral perturbation theory. This led the authors in ref. [3] to perform fits to Nf = 2 chiral 
expansions. The results obtained in ref. [3] that are relevant for our discussion are summarised in 
table El 

A latter article of the same collaboration |54j considered simulations directly performed at the 
physical point including non-perturbative renormalisation. This has induced a significant modifica- 
tion for the quark mass renormalisation factor, becoming Z m = 1.441(15) (non-perturbative) instead 
of Z m = 1.114 (one-loop perturbation theory) leading to an increase (decrease) in the values of quark 
masses (condensates) by a factor 1.30. This should be taken into account when comparing the results 
obtained from the PACS and RBC/UKQCD sets in this article. The results obtained in ref. [54] that 
are relevant for our discussion are summarised in table 01 Since the simulation was performed at the 
physical point, there is no further information on LECs describing the pattern of Nf = 2 and Nf = 3 
chiral symmetry breakings. 

r I 31.2 ±2.7 

F K /F n 1.333 ±0.072 

m s (2 GeV)[MeV] 92.75 ± 0.58 ± 0.95 
m(2 GeV) [MeV] 2.97 ± 0.28 ± 0.03 



Table 4: Results obtained by the CP-PACS collaboration [54] with non-perturbative renormalisation 
and simulation at the physical point. 



B Dimensional estimate of the remainders 

As indicated in eqs. (185j) -(|87j). we can estimate the higher-order remainders assuming that they are 
dominated by NNLO contributions and using resonance saturation. A typical order of magnitude for 
0(m1) remainders so that the chiral series converge is 10% so that Ah — 0.8 GeV. The corresponding 
size a of the remainders is given in table 02 and the remainders will be required to stay in the range 
[—a, a] in our fits to lattice data if necessary. In the specific case of the electromagnetic square radius 
of the pion, we have combined the uncertainty on the experimental measurement of the square radius 
with the theory uncertainty on the remainder in quadrature (the range for the kaon radii would be 
the same). 



Remainder 


G 


d, e, d + 


0.148 


d',e',d- 


0.024 


e+ 


0.005 


e_ 


0.001 


e v 

°7T 


0.318 



Table 5: Size of the NNLO remainders allowed in our fits, based on a dimensional estimate. 
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C Fits with "Total" data sets 

In Sec, 16.21 we considered fits to both PACS-CS and RBC/UKQCD data restricted to the low-mass 
and low- momentum region ("Subset" data). We have also performed fits to the whole sets of data 
available ("Total" data), in order to test the stability of our results, and to illustrate the interest 
of having larger data sets to determine NNLO remainders in an accurate way. We are aware that 
some of the data points considered here may stand outside the region of validity for Re^PT, but 
we found nevertheless interesting to provide these results, showing a good consistency with those 
obtained with "Subset" data. 

Our results are summarised in Tables [6] and [71 The first series of rows corresponds to the outcome 
of the fit, whereas the lower rows are quantities derived from the results of the fit (LO LECs, NLO 
LECs, quantities in the Nj = 2 chiral limit, quantities, relative fraction of LO/NLO/remainders 
contributions at the minimum for several observables), and the last row is the x 2 P er degree of 
freedom. Most of the comments made in Sec. 16.21 can be restated, with a few changes in the case of 
the RBC/UKQCD data (larger value of Fx/F n and lower value of /+(0) than in the "Subset" case). 
We notice that the fits are fairly good, and that all NNLO remainders turn out to lie within their 
expected range. 
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0.88 ± 0.01 


0.87 ± 0.01 


5(2 GeV)[GeV] 


2.34 ± 0.07 


3.58 ± 0.11 


F[MeV] 


86.7 ±0.3 


86.3 ±0.3 




3.1 ±0.6 


3.8 ±0.6 




3.9 ±0.2 


4.2 ±0.2 




2.07 ±0.15 


1.88 ±0.14 


B/B 


1.32 ± 0.10 


1.35 ± 0.11 


F/Fq 


1.25 ± 0.04 


1.18 ± 0.04 


/+(0) 


1.006 ±0.149 


1.011 ±0.149 




0.56 ±0.54- 0.10 


0.63 ±0.28 ±0.09 


*1 


0.38 ±0.69 -0.07 


0.41 ±0.52 ±0.07 




0.44 ±0.67- 0.11 


0.48 ±0.50 ±0.02 




0.31 ±0.77- 0.08 


0.33 ±0.65 ±0.02 


x 2 /iv 


13.6/7 


13.8/15 



Table 6: Results of fits performed on the data from the RBC/UKQCD [S EJ [6] and PACS-CS [3] 
collaborations on pseudoscalar masses and decay constants, considering all the available unquenched 
data (Total). Only statistical errors are shown. The LECs are given at the scale (i = m p . The 
-fQ3 form factor at zero momentum transfer is a prediction of the fit (with an error combining those 
obtained from the fit and the maximal contribution allowed for the remainder from dimensional 
estimation). The penultimate set of rows collects the relative fractions of LO/NLO/remainders for 
decay constants and masses at the minimum. 

[16] B. Ananthanarayan, G. Colangelo, J. Gasser and H. Leutwyler, Phys. Rept. 353, 207 (2001) 
[arXiv:hep-ph/00052"97] . 



28 



With K n 
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RBC/UKQCD Total II 
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m s {2 GeV)[MeV] 
m(2 GeV)[MeV] 
S (2 GeV)[GeV] 
F [MeV] 
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Table 7: Results of two different fits of the data from the RBC/UKQCD [2J [5l [6] on pseudoscalar 
masses and decay constants, as well as on form factors. We considered all the available un- 
quenched data (Total), and either fixed the lattice strange quark mass (marked then with a star) or 
let it vary freely. Only statistical errors are shown and LECs are given at the scale \i = m„. The 
penultimate set of rows collects the relative fractions of LO/NLO/remainders at the minimum for 
decay constants, masses and form factor at vanishing transfer momentum. 
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